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rj^ We review the current understanding on universal behaviors in granular flows through 

a vertical pipe and traffic flows. We carry out weakly nonlinear analysis of a model for 
traffic flows based on the technique of soliton perturbations, and determine the selected 
propagating velocity, the amplitude, the width of interfaces connecting between jam phase 
and non-jam phase. From the direct simulation of the model, we have confirmed the validity 
of our theoretical analysis. We also introduce a model for granular pipe flow supplemented 
by the white noise, which reproduces P(f) ~ /~ 4 ^ 3 , where P(f) is the power spectrum in the 
^ , frequency /. (This paper will be published in Progress of Thoeretical Physics Supplment). 

7— I . 

^ §1. Introduction 

Q\ , Recently, cooperative dynamics in dissipative systems consisting of discrete ele- 

ments have attracted much attention. Researches on granular materials are efforts 
to understand unusual behaviors of discrete element systems© such as convection!) , 
size segregation EP, bubblinglf, standing waves & and localized excitations as well 
as thermodynamic descriptions of granular particles under the vertical vibrationsEP . 
In particular, it is interesting that jam formation of particles and a power law in the 
q \ power spectra in flows of granular particles through a narrow vertical pipe. Simi- 

larly, traffic jams in a highway is also an attractive subject not only for engineers 
but for physicists Similarities between two phenomena are obvious. Both consists 
of discrete dissipative elements, vehicles and particles which are confined in a quasi 
one-dimensional systems such as a highway and a pipe. There is an optimal velocity 
in each system, the competition between the relaxation to the optimal velocity and 
acceleration of particles produces jam formation. We, thus, expect that there exists 
common and universal mathematical structure behind these phenomena. 

In the next section, we will introduce typical models for granular flows and traffic 
flows, and summarize the current status of our understanding on universal properties 
in one-dimensional flows. Among these studies we will focus on recent two main 
streams in studies of granular flows and traffic flows. First, we will review the recent 
progress in theoretical analysis of pure one-dimensional models. The second, we 
will introduce the progress in analysis for power spectra in density auto-correlation 
function in quasi one-dimensional systems. 
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From the theoretical analysis of pure one-dimensional models, at least, it has 
been confirmed that perturbed solitons by dissipative corrections play important 
roles in traffic flows and granular flows. In particular, Komatsu and SasaB) have re- 
vealed the mechanism of jam formation iri a traffic flow by the perturbative treatment 
of solitona. Hayakawa and Nakanishill2P have generalized the analysis of Komatsu 
and SasacP, and demonstrate that universal mathematical structure exists in pure 
one-dimensional models for granular flows and traffic flows. In section 3 consisting 
of three subsections, we will review the details of theoretical argument by Hayakawa 
and Nakanishili2P which discussed a model of traffic flow. We will stress that the 
framework of our analysis can be used in any one-dimensional models for granular 
flows and traffic flows. 

In realistic situations, however, e.g. highways have several lanes, and vehicles 
(particles) can pass slow vehicles (particles). When we include multi-lane effects 
in one-dimensional models, separations between jam and non-jam phases become 
obscured. However, there is another universal law in quasi one-dimensional systems 
for flows of dissipative discrete elements, i.e. a power law in power spectrum of 
density auto-correlation function of vehicles or particles. Recently, Moriyama et 
al.liiP have presented that the power spectrum is given by P(f) ~ f~ a with a = 
1.33 from their experiment on granular flow through a vertical pipe. This result is 
expected to be a universal in quasi one-dimensional dissipative flows such as traffic 
flows in a highway. Thus, we will clarify the mechanism to appear /" 4 / 3 law m 
section 4. For this purpose, we will introduce a simple model supplemented by the 
white noise. From the simulation of our model we will confirm that our model can 
reproduce / _4//3 . We will explain the mechanism how to obtain /~ 4 / 3 -law from the 
simple analytic calculation. 

In section 5, we will give concluding remarks, especially on universality in gran- 
ular flows through a pipe and traffic flows. We also summarize our results. 

§2. Models 

In this section let us summarize what models exist and what consensus in studies 
of granular flows and traffic flows is obtained. 

There are many models to describe traffic flows and granular flows through a 
pipe. We believe that universal behaviors do not depend on the choice of the model. 
Recently, Hayakawa and NakanishiEf have proposed a generalized optimal velocity 
model for traffic flows 

x n = a[U(x n+ i - x n )V{x n - x n -{) - x n ], (2-1) 

where x n and a are the positions of n th car, and the drag coefficient, respectively. 
This model contains the psychological effect of drivers. Namely, the driver of x n takes 
care of not only the distance ahead x n+ \ —x n but also the distance behind x n — x n -\. 
The optimal velocity function U should be a monotonic increasing function of the 
distance of x n+ \ — x n and V should be a monotonic decreasing function of x n — x n _\. 
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Thus, we adopt 

U(h) =tanh(/i-2)+tanh(2); V(h) = 1 + / (1 - tanh(/i - 2)) (2-2) 

for the explicit calculation in section 3, where /o is a constant. We put these optimal 
velocity functions as the product form UV in ( |2-1| ), because the driver of x n cannot 
accelerate the car without enough space ahead even when the distance x n — x n —i 
becomes short. In other words, the model including U(x n+ i — x n ) + V(x n — x n -i) is 
mathematically unstable and unphysical, because the acceleration by V causes crash 
of vehicles. This model (2-1) with (2-2) is a generalization of the optimal velocity 



(OV) model proposed by Bando et al. tn> 

x n = a[U(x n+1 - x n ) - x n \. (2-3) 

The generalized OV model is similar to the model of granular flow in a one- 
dimensional tube 



C[U(x n+ i - x n -i) - x n ] + T[ip'(x n+ i - x n ) - ip'{x n - X n -l)], (2-4) 



where £ and T are respectively the drag coefficient and the strength of collision 
among particles. U and ip are the optimal velocity and soft core repulsion potential, 
respectively©. 

There is a fluid field model to describe traffic flowst^P which consists of mass 
conservation for density field p and momentum conservation for velocity field v as 

d t p = -d x (pv), 

dt v = -vd x v - ^d xP + Up{p) ~ V + 2^ v> (2-5) 
P t p P 

where t p and T e are a characteristic time for the relaxation and the effective temper- 
ature, respectively. This model also contains the relaxation mechanism to an optimal 
velocity U p , while the pressure term and the viscous term are phenomenologically 
introduced to stabilize the solution of a set of equations (2-5). 

This fluid model (2-5) i^imilartoiluid models to describe granular flows through 
a pipe and fluidized beds E-§' iHir t3)'EZP and mixture of polymersll^P. As a typical 
example, we write an explicit form of a fluid model for granular flow by Sasa and 
Hayakawa© at the Froude number Fr: 

d t p = -d x (pv), 

8 t v = -vd x v-((p)(V -v)--± r -f"(p)d x p + Kd 3 xP 

b r 

+-d x (pu{p)d x v). (2-6) 
P 

Although this model seems to be complicated, each term corresponds to that in 
eq.(2-5). Here the term in proportion to k is introduced to stabilize the solution 
furthermore, which can be regarded as a coupling term in mixing free energy. The 
pressure f'(p) arises from the collisions among particles, which is a recovery force to 



4 



the mean density po and diverges at the closest packing states at p cp . The kinetic 
viscosity v (/?) also diverges at p cp . Komatsu and Hayakawa^, thus, adopted f(p) = 
(3\{p — po) 2 /(p C p — p) and v(p) = fo/iPcp — p) for the simulation of (2-6), where (5\ 
and are constants. The optimal velocity Vq is assumed to be a constant because 
of the incompressibility of the fluid, and is a phenomenological function to be 
determined by the sedimentation rate. 

At first sight fluid models (2-5) and (2-6) are very different from discrete models 
such as ( p-ip , ( [2 -3D and (2-4). However, there is common mathematical structure. 
The fluid models of granular flows such as eq.(2-6) are reduced to the Korteweg-de 
Vries (KdYl equation near the neutral curve of the linear stability©' §)•§). Kurtze 
and HpngliH) also derived KdV equation from the fluid model of the traffic flow 
(2-5)113). Of course, it is easy to derive KdV equation from the discrete models 



(2_1), (2^3) and ( p-4[) near the neutral curve. Thus, at least, we get a consensus that 
dissipative solitons play important roles for granular flows and traffic flows. 

Unfortunately KdV equation is not adequate to describe the traffic jams, be- 
cause solutions of these equations are essentially pulses and no interface solutions 
connecting jam phase with non-jam phase are included. Komatsu and S solved 
such a puzzle from the analysis of the original OV model Q2-3| ). They have showed 
that fl2-3| ) can be reduced to the modified KdV (MKdV) equation at the critical 
point (the averaged car distance h = 2) or the most unstable point on the neu- 
tral curve. They also show that symmetric kink solitons deformed by dissipative 
corrections describe a phase separation between bistable phases. This analysis is 
also consistent with recent analysis for the exactly solvable models which may be 
regarded as simplified optimal velocity models®. However, as will be shown, the 
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generalized optimal velocitv model (^J) and granular model fl2-4| ) as well asthe 
model (2-5) of traffic flowst$ and fluid models, e.g. (2-6) for granular flowst^P't 
are not reduced to MKdV equation at thecritical point or the most unstable point 
on the neutral curve. In fact, KomatsutlP has shown that the fluid model (2-5) 
has the following properties: (i) Interfaces (kinks) between jam and non-jam phases 
are asymmetric, (ii) The critical point to appear kinks is, in general, different from 
the most unstable point on the neutral curve, (iii) Eventually one branch of the 
coexistence curve exists in the linearly unstable region. He also demonstrates that 
MKdV equation is recovered in a special choice of parameters of the fluid model, 
while fluid models cannot be reduced to MKdV equation in general cases. Thus, we 
need to clarify universal characteristics of dissipative particle dynamics in general 
cases which contains (2-1), (2-4) and fluid models i3' v3' 119' i$ . For this purpose, 
we will focus on analysis the simplest model (|2-1| ) among them to characterize the 
phase separation between jam and non-jam phases. 

We may be suspicious of reality and relevancy of pure one-dimensional models. 
As indicated in Introduction, when we include multi-lane effects in one-dimensional 
models, separations between jam and non-jam phases become obscured. However, 
the fundamental characteristics of pure one-dimensional models should be impor- 
tant for quasi one-dimensional cases. To demonstrate the relevancy of pure one- 
dimensional models, we will focus on another universal law in quasi one-dimensional 
systems, i.e. a power law in power spectrum of density auto-correlation function of 
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vehicles or particles. 

Power-law form of the power spectrum P(f) ~ f a , whepe / is frequency, of 

tnsity fluctuations was found in both numerical simulations!^ and experiments 
''EiP. Although their interpretations on the origin of the emergence of density 
waves are different, estimated values aLthe exponent a is close to each other (1.3 < 
a < 1.5). The previous reports E^'EJ on the estimation a = 1.5 seem a little 
ambiguous: The volume of air flow out of the bottom end of the pipe was not well 
controlled. Besides, the power spectra they obtained were still noisy. Recently, 
Moriyama et al.lHP have presented better-controlled air flow out of the pipe and 
more accurate experimental results than the previous ones by increasing the number 
of trials. One of their results is the precise estimation of the scaling exponent of 
the piwer spectrum P(f) ~ f~ a with a = 1.33. This result is identical to that by 
LGAcHP, and is expected to be a universal law in quasi one-dimensional dissipative 
flows such as traffic flows in a highway. 

The second purpose of this paper is to clarify the mechanism to appear /~ 4 / 3 law 



in power spectra. We, thus, extend the one-dimensional model (2-4) to a stochastic 
model supplemented by the white noise as 

K = C[U( hn+ \ +hn ) - U{ *±±±=L )) - h n ] 

+T[ip'{h n+1 ) + <p'(hn-i) - 2ip'(h n )] + f n (t), (2-7) 

where h n = x n +i — x n . The most crucial simplification of the model is that f n is 
assumed to be the Gaussian white noise with zero mean. The optimal velocity U(h) 
may be the sedimentation rate which is a function of the local volume fractiont^ 
in general. It should be noticed that the drag £ is irrelevant in systems where the 
bottom end of the pipe is fully open, because air in the pipe flows away together 
with particles. Thus, to observe density waves it is important to close the cock of 
the pipe. 

One of the important points is that the white noise is introduced not to an 
equation for x n but to an equation for the relative motion of particles h n , because 
the ordering of particles along a pipe or a highway is not conserved due to passing 
of particles in multi-lane or multi-dimensional systems. From both simulation and 
analytic calculation of (2-7), we will demonstrate that this simple model reproduces 
a = 4/3 near the neutral curve of the linear stability analysis of uniform states in 
section 3. We can expect the existence of a universal law regardless to the choice of 



a specific model. Therefore, the model (2-1) for h n supplemented by the white noise 



is expected to belong to the same universality class that by (2-7). The fluid models 
(2-5) and (2-6) with non-conserved white noise also should behave similarly as those 
for discrete models. 

Now, we have come back to an interesting question: What are good models 
among many models? The answer is simple. The simple models are good ones. 
For example, when we compare discrete models ( pT[ ) and (2-4) with fluid models 



(2-5) and (2-6), the representation of discrete models are shorter than fluid models. 
Since fluid models are partial differential equations, we may need careful check of 
the validity of the simulation scheme and long CPU time to simulate them. On the 
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other hand, discrete models which are coupled ordinary differential equations are 
free from such problems. The theoretical analysis for discrete models is also simpler 
than that for fluid models, because each discrete model contains only one set of 
equations. On the other hand, eq.(|2-3j) is a oversimplified model which loses some of 



universal properties. Therefore, we believe that models (2T) and fl2-4| ) are the most 
fundamental ones. In the following sections, thus, we will focus on the analysis of 
discrete models fl2-l|) and (2-7). 



§3. Theory of a pure one dimensional model for traffic flow 



In this section, we concentrate on the analysis of ( |2-1| ). The result explained in 
this section is universal for all of models introduced in the previous section. This 
section consists of three subsections. In the first subsection, we will summarize the 
result of linear stability analysis for uniform flows. In the next subsection, the main 
part of this section, the details of weakly nonlinear analysis will be explained, where a 
steady propagating solution is presented and the selection of its propagating velocity, 
the width of interfaces and the amplitude will be discussed. In the last subsection, 
we will confirm the validity of our theoretical analysis from the comparison between 
direct simulation and theoretical results. 

3.1. Linear stability of uniform flow 

In this subsection we summarize the linear stability analysis of the uniform 
propagating flow. 

It is obvious that there is a constant propagating solution with x n+ \ — x n = 
constant. Let us rewrite (|2-1| ) as 

r n = a[U(h + r n+1 )V(h + r n ) - U{h + r n )V{h + r„_i) - r„], (3-1) 

where h is the averaged distance of successive cars and r n is x n+ \ — x n — h . The 
linearized equation of (|3-1|) around r n (t) = is given by 



r n = a[U'(h)V(h)(r n+1 - r n ) + U{h)V(h)(r n - r„_i) - f n ], (3-2) 

where the prime refers to the differentiation with respect to the argument. With the 
aid of the Fourier transform 

1 N 

r q (t) = — exp[-iqnh]r n (t) (3-3) 
n=l 



with q = 2irm/Nh and the total number of cars iV we can rewrite ( |3-2[ ) as 

(d t -a + (q))(d t -a^q))r q (t)=0 (3-4) 

with 

a ± (q) = -| ± VW2) 2 - aD h [U,V](l - cos(qh)) + ia(UV)' sm(qh) , (3-5) 

where we drop the argument h in U and V. Dh[U,V] = U'(h)V(h) — U(h)V'(h) 
denotes Hirota's derivative. The solution of the initial value problem in Q3-4j) is the 
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linear combination of terms in proportion to exp[o"+(g)i] and exp[cr_(g)t]. The mode 
in proportion to exp[<r_((/)i] can be interpreted as the fast decaying mode, while the 
term in proportion to exp[cr+((/)t] is the slow and more important mode. 

The violation of the linear stability of the uniform solution in ( |3-2| ) is equivalent 
to Re[a+(q)] > 0. Assuming qh ^ (qh = is the neutral mode), the instability 
condition is given by 

2(in/) /2 cos 2 (^) >aD h [U,V]. (3-6) 

Thus, the most unstable mode exists at qh — > and the neutral curve for long wave 
instability is given by 

2(UV)' 2 

The neutral curve in the parameter space (a, h) is shown in Fig.l for /q = 1/(1 + 
tanh(2)) in fl2-2j ). For later convenience, we write the explicit form of the long wave 
expansion of <r+ in the vicinity of the neutral curve 

a + (q) = lCo qh - ^°-^{qhf - l ^fc, - ^-c 2 + 0((qh)% (3-8) 
a n [hy 6 4:a n (h) 

where Co = (UV)' . Thus, the uniform state becomes unstable due to the negative 
diffusion constant for a < a n (h). 

3.2. Nonlinear analysis 

The simplest way to describe nonlinear dynamics is the long wave expansion 
with the help of a suitable scaling ansatz. It is easy to derive the KdV equ ation 
near the neutral curve from (2-1) as in the case of fluid modelsli^'Ej^EZP'll^. As 



mentioned in Introduction, to describe the phase separations, however, we should 
choose the critical point (a,h) = (a c ,h c ) at (U(h)V(h))" = where the coefficient 
of d x r 2 becomes zero on the neutral curve. At this point the cubic nonlinear terms 
can produce the interface solution to connect two separated domains. The explicit 
critical point of (|2-2j ) with /o = 1/(1 + tanh(2)) is given by 



512 

he = 2-tanh _1 (l/3) ~ 1.65343; a c = — / 2 ~ 1.63866. (3-9) 

81 

Unfortunately, the reduced equation based on the long wave expansion of our model 
is an ill-posed equation. In fact, the scaling of variables as r n (t) = er(z,r), z = 
e(n + cot) and r = e 3 i with e = \J (a c — a)/a c leads to 

8 T r = ai d z r 3 - a 2 d 3 z r + a 3 d 2 r 2 (3-10) 

in the lowest order, where a\, a 2 and 03 are constants. The solution of ( |3-10| ) is blown 
up within finite time. The reason is simple. Its linearized equation around r = d$ 
is unstable for all scale, because the solution with r — do ~ exp[ikz + A^r] has the 
growth rate i?e[Afe] = 2k 2 a^alQ which is always positive when m^oIq > 0. Therefore, 
the simple long wave expansion adopted by Komatsu and SasaB-* for (|2-3| ) and 03 = 
cannot be used in our case. 
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This short scale instability in ( [3-10 ) arises from the long wave expansion. To 
avoid such the difficulty, we only focus on a steady propagating solution for the 
theoretical analysis. Dynamical behavior is followed by the direct simulation of 
discrete model ( |2T| ) , where simulation of ( |2-1[ ) is much easier than that of the reduced 
partial differential equations such as MKdV. In addition, our long time simulation 
suggests that the solution of (|2-l| ) seems to be relaxed to a steady propagating mode. 

To obtain the steady propagating kink solution, at first we eliminate the fast 
decaying mode in (|3-1|) as 



(d t -a + (d x ))r(x,t) = (a + -a^)- 1 N[r(x,t)}, (3-11) 
where N[r] represents the nonlinear terms 

= U(h + e hd *r)V(h + r) - U(h + r)V(h + e~ hdx r) 

a 

-U'(h)V{h)(e hd * - l)r(x,t) + U{h)V\h){\ - e - h9 *)r(x,t).(3-12) 



Since (<r+ — cr_) 1 is the inverse of the polynomial of the differential operators, it is 

convenient to use the expansion (a + —a)~ 1 ~ a _1 [l (UV)' d x +0(h 2 )]. It should 

i 1 a 

be noticed that ( |3-11| ) contains most of important information and no instability in 

time evolution. 

To obtain the scaled propagating kink solution we assume the scaling of the 
variables by e = \J [a c — a)/a c as 



r(x, t) = ^ )ml R(z) : z = e^{n + c t - e 2 j(t)t), (3-13) 

where the arguments in U and V are fixed at h = h c , and 7 is the positive free 
parameter which will be determined from the perturbation analysis. The expansion 
of N[r] is given by 

00 00 

N[r]/a = E E h m C mn d™r m - h 3 U'Vd x rd 2 x r + • • • , (3-14) 

n=l m=2 

where C 21 = \{UV)" ,C 22 = \D h [U, V}' , C 23 = ±{UV)",C 3 i = \{UV)"> \C 32 = 

±D h [U,V]",C 41 = ±(UV)"" with D h [U,V]' = -^D h [U,V}. Substituting O 
into ( 3- 11| ) with the help of ( |3T4 ) we obtain 

t It? - Ri ~ R2 ~ + T^) = e T M ^ ^ 

dz I dz A dz I dz 

where = 3D h [U, V]' '/ '(2y 'cq\{UV)"'\) and 

,, rn1 ^ (dR\ 2 dR 3 ni 1 t .dR d 3 R 2dR„ . r zR 

(3-16) 
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Here 7 = dj/dr with r = et, and I/77 = y/6D h [U, V]/c , p 23 = 3V6U'V' / ^c \(UVy 
P32 = V3/2D h [U,V}"/\(UV)'"\ and p u = y/3^(UV)""/(2y/2\(UV)>"\*). 
Assuming R(z) = Rq(z) + eRi(z) + • • we obtain a solution 



R^\z) = tanh(9±z); 0± = - - ^- + 2 (3-17) 

in the lowest order. This solution represents a kink or an anti-kink connecting 
between jam and non-jam phases because of 9 + > and #_ < 0. Notice that the 
solution is not localized one and does not satisfy the periodic boundary condition. 
Therefore, we neei a careful treatment of the boundary condition. In fact, our 
preliminary resulted suggests that the selected value of 7 under the open boundary 
condition is different from that under the periodic boundary condition. 

In this paper we restrict ourselves to the case under the periodic boundary. To 
satisfy the periodic boundary condition, we use 

R (z) ~ R {+) (z - z+ )-l + Ro { ~Hz - z-) (3-18) 

as an approximate solution of the lowest order equation ( |3-15| ) , where a pair of kink 
and anti-kink exists at z = z + and z = z~. Since ( 3T8| ) is not an exact solution of 
( 3T5| ), there should be an interaction between the kink and the anti-kink which is 
an exponential function of the distance between themt^P. 

Now, let us discuss the effect of perturbative terms in ( |3T5 ). It is known that 
perturbation of solitons or solution including a free parameter b ec omes unstable ex - 
cept for the solution where the parameter has a special value. B^'EZP'E3'Ej)'Ej)'E1P'E3. 
The linearized equation of ( |3-15 ) can be reduced to 

CR 1 = ^-M[R ], (3-19) 

dz 

where 

C = &l + d z - 6R d z - 3R%d z + 2pd 2 z R + 4(3d z R d z + 2[3R d 2 z • (3-20) 

To obtain a regular behavior of perturbation in O(e) the perturbed solution should 
satisfy the solvability condition 

(!Pb, ^M[Ro}) = lim [ L dz%-^-M[R ] = 0, (3-21) 

dz L^qoJ^l dz 

where L is the system size and •f'o satisfies 

tf& = 0: C ] = -d 3 z -d z + 3R 2 l d z + 2(3R d 2 z . (3-22) 

When we adopt ( |3T8|) as Rq, &o also should satisfies the periodic boundary condition. 
Thus, we assume 

%(z) = V {+) (z - z+) - 1 + ^-\z - *_), (3-23) 

where is the solution of ( 3-22|) when we replace Ro by Rq ■ Since ^0 = d z &o 
satisfies 

&$o(z) =0; C) = -d 2 z - 1 + 3R\ + 2(3R 8 Z , (3-24) 
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the solution of ( 3-22| ) can be expressed by 



where we use 



a± 

~2 



dz / (sech[0±* / ]) 1 /*±; ^ ±] {z) 



-iPi 



(±), 



(3-25) 



4%) = (sech[6>±z]) 1/y ± 



V«3 



The constant a>± in ( |3-25| ) is determined to satisfy \Pq(±oo) 



(3-26) 
-1. Thus, we obtain 



a± 



20± 



j(±) 

1 r> 



dx(sechx) 



i/ei+2n 



r(l/(2f£)+n) 
r(l/(20|)+n + l/2) ! 



(3-27) 



where -T(x) is the gamma function. It should be noticed that ^q{z) is n °t a local- 
ized function. Therefore, we cannot neglect the boundary effects in the solvability 
condition (3-21). 

Let us rewrite ( f3-21| ) as 



($ (z),M[Ro\), 



(3-28) 



where [f(z)]~L = f(L) — f(—L). From (3T6) it is obvious that contribution from 
terms except for those in proportion to p4i and zR is zero in the left hand side of 
( 3-28; ) under any boundary conditions. Notice that the contribution from the term 
in proportion to tR vanishes because of its symmetry. If we adopt the periodic 
boundarycondition and use (3-18) and fl3-23|) , the contribution from the term is 
canceledH. Thus, the left hand side of (|3-28j) is reduced to 



[%M[R ] 
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5/2 



L. 



(3-29) 



On the other hand, the right hand side of ( |3-28| ) is the integration of the product of 
d34(i| ) and ( gjjj ). 

For simplicity, to obtain the explicit form we assume fo = 1/(1 + tanh(2)) in 
(|2-2|). In this case the coefficients in (3-31) are reduced to P23 = —3/2, P32 = —f3, 
p 41 = -1/4, T) = 1/(4/3) c = 2 6 /o/3 3 = 1.20689, and (3 = 3\/3/(8\/2/o) = 0.902037, 
6co/\(UV)'"\ = 9/4, 8+ = 1.2897187, 6L = -0.387814. Finally we obtain 



{L-(9 + -6_)}^ = Af3 1 2 { w 
(±) 



+ 1 



7+ 



+ 1 



7- 



(3-30) 



where we use 



I 



n+l 

r(±) 
in 



2n6 z ± + 1 
(2n + 1)9 2 ± + 1 



. Here j± is given by 
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Thus, we obtain j± as 

_ 8f36±(36 2 ± + 1) 
7± ~ 37^ + 8^-8" 
Substituting ( ggg) into (|3~30| ) we obtain 



7 



+ e_)^ = ar «(l_-L) j 



(3-32) 



(3-33) 



where 



a 



24/3 yffiT^ 
4/3 2 + 9 ' 



7 



3(12/? 2 + 25) 
61/3 2 + 132 



0.574189' 



(3-34) 



From ( 3-13 ) the amplitude Ae of Rq can be regarded as the order parameter of phase 
separation, which is given by 



= 2 e V7* 



1.13663e. 



In the vicinity of 7* the time evolution of 7 is described by 

07* 

7(7-) ~ 7* + Aexp[ — t] 

1j 



(3-35) 



(3-36) 



Notice that 7* is the stable fixed point in the time evolution ( |3-33| ). where we use 

Two remarks on the result of this section are addressed: We recall that to 
derive ( |3-34 ) from ( 3-28| ) we assume /o = 1/(1 + tanh(2)). Although to obtain the 
result for any /o is not difficult, we omit such a generalization to avoid long and 
tedious calculations. Thus, the expression with [3 = in (g]34|) does not recover the 
result by Komatsu and SasaEP. As the second remark, the characteristic time for 
the relaxation ( |3-33| ) or ( |3-36| ) is proportional to the system size L. This result is 
reasonable, because the time needed to simulate (|2T|) is proportional to number of 



cars N. This tendency has been confirmed by our simulation. 



3.3. Simulation 



To check the validity of our analysis in the previous subsection we perform the 

id fl2% with f Q = 1/(1 + tanh(2)) near the critical 
under the periodic boundary condition. We adopt the classical fourth- 
order Runge-Kutta scheme with fixed time interval At = 2~ 4 . Since our purpose is 



numerical simulation of (2T 
point 



the quantitative test of (3-17) and ( |3-34| ), the initial condition is restricted to the 
localized symmetric form r n = 18.7/A r (tanh(n — N/4) — tanh(n — 3A r /4) — 1) where N 
is the number of cars. Taking into account the scaling properties we perform the sim- 
ulation for the set of parameters (e,N) = (1/2, 32), (1/4, 64), (1/8, 128), (1/16, 256) 
until r n relaxes to a steady propagating state. Our results are plotted in Figs.l and 
2. 

Figure 1 displays points which have the maximum h max and the minimum h m i n 
values of successive car distance in each parameter set, and theoretical coexistence 
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curve 



1 



(h - hcf 
A 2 



A = 1.13663 • 



(3-37) 



where we use a = a c (l — e 2 ) and Ae = h — h c . The agreement with each other 
is obvious. From this figure we can see that one of the branches is in the linearly 
unstable region but the theoretical curve recovers the simulation result. We stress 
that the evaluated value of A from simulation at e = 1/16 is 1.14273 • • •. Thus, the 
deviation between simulation and theory is only 0.53 %. 

1.8 
1.7 
1.6 
1.5 



1.4 
1.3 
1.2 
1.1 







i i i i 

co(h) 






an(h) 




/ y 


'data' » 




o I 





0.6 0.8 1 1.2 



1.4 1.6 
h 



1.8 



2.2 2.4 



Fig. 1. Theoretical coexistence curve (solid line) a = a c (1 - {h - h c f /A 2 ) where A = 1.13663, 
h c = 1.65343 and a c — 1.63866, and the neutral curve (broken line) for the model described 
by eq.(2) with W(h) = tanh(ft - 2) + tanh(2), V(h) = 1 + (1 - tanh(ft - 2))/(l + tanh(2)) 
and f n — 0. h denotes the average distance between successive cars. The data is obtained for 
minimum and maximum values of r n at a given a. 10 ^ 



Figure 2 demonstrates that the numerical result has a scaling solution which has 
an asymmetric kink-antikink pair. The linear combination of our theoretical curve 
( 3- 17] ) with ( 3-32 ) is plotted as the solid line by choosing the position of the kink and 
the anti-kink. Our theoretical curve agrees with the result of our simulation without 
other fitting parameters. Thus, we have confirmed the validity of our theoretical 
analysis. 



§4. / 4 / 3 law in power spectra 



The purpose of this section is to clarify the mechanism to appear P(f) ~ / 4 ^ 3 
law in power spectra. We, thus, extend the one-dimensional model ( |2-4| ) to a stochas- 
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-1.5 1 1 1 1 1 u 

50 100 150 200 250 

z 

Fig. 2. Theoretical curve (solid line) and scaled data obtained from our simulation for scaled r n . 
Each data denotes (e, N) = (1/2, 32), (1/4, 64), (1/8, 128), (1/16, 256), where e = (1 - a/a c ) 1/2 , 
'N.s' represents the data for N cars (particles). Theoretical curve is given by R(z) — 
tanh(£0 + (z-z + ))-l+tanh(£0_(z-z_)) with £ = (67*) 1/2 /16, 7 * = 0.574189 , 9+ = 1.2897187 
and 9- = —0.3876814, where only z + — 62.5 and Z- = 190.5 are fitting parameters. Spatial 
scale is measured by the scale for N = 256. ^ 



tic model (2-7) supplemented by the white noise. We will demonstrate the simple 
model reproduces a = 4/3 near the neutral curve of the linear stability analysis of 
uniform states. 

At first, we briefly summarize the result of linear stability. Although the frame- 
work is common with that in section 3.1, the explicit forms are a little different. The 
neutral curve is given by T n = U' 2 /<p"(h). The critical point is given by the cross 
point of If" (h c ) = on the neutral curve. We will focus on the behaviors for weakly 
unstable or stable region at T = T c (l — [/,) with <C 1. 

Let us simulate (2-7) directly Adopting 

U(r) = tanh(r - 2) + tanh(2), ip(r) = sech 2 (r) (4-1) 

with £ = 2, iV = 256, T c = 3.95798 • • • , and h = 2 at t = 0, we numerically 
calculate (2-7) by the classical Runge-Kutta method until t = 2 11 with time interval 
At = 1/2 4 under the periodic boundary condition. We use the uniform random 
number distributed between —X and X with X = 9/1024 for f n (t). Figure 3 displays 
the power spectrum P(f) =< \p(f)\ 2 > obtained from our simulation of (2-7) at \i = 
1/64, where p(f) is the Fourier transform of the discretely sampled data of the density 

p(t) = — ^ — -— with the time interval 1. This clearly supports P(f) ~ /~ 4 ^ 3 law 
A n r n (t) 
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in the range of / between 10 and 1000 as in the experimenting. It should be noticed 
that the data in Fig. 3 may suggest steeper slope than /~ 4 / 3 in low frequency region 
/ < 10. Although the steep slope close to / _3//2 at small / is not observed in the 
experiment by Moriyama et al. t3 , it can be explained easily by the diffusive behavior 
of an interactive pair of kink and antikink. The details of the process to produce 
/ -3 / 2 will be reported elsewhere. 

From the examinations of several values of fj,, we have confirmed that the quali- 
tative results are insensitive to the sign of p when \p\ <C 1. This result is reasonable 
because near the neutral curve the time scale of relaxation or growth of fluctuations 
is much longer than the time scale induced by the noise f n {t)- Our numerical re- 
sult suggests that the linear relaxation theory of fluctuations can be used to explain 
P(f) ~ /~ 4 / 3 - This / _4//3 law can been seen from the simulation of traffic flow (2-1 
with the noise. 



100 



10 



0.1 



P(f) 

'DATA' 



0.01 



10 



100 



1000 



Fig. 3. Log-log plot of power spectrum P(f) obtained from the simulation of eq.(2.7) with (4.1), 
where the unit of / is l/(27r) and the unit of P(f) is not normalized. The solid line represents 



Thus, let us briefly explain how to appear / 4 / 3 law from the behavior of struc- 
ture factor 



Sk(t) =J2< exp{ik{r n (t) - r m (0))] >= — ^ exp 



2 0rwn(v 



(4-2) 



where nTO (i) =< (r n (t) — r m (0)) 2 >. Notice that the structure factor is directly 
related to the auto-correlation function as Sk(t) = Jq dt' < pk(t + t')p^{t) >, where 
Pk(t) is the Fourier component of the density field of particles p(r, t). In weakly stable 
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states, i.e-—\x < and <C 1, Sk(t) can be calculated as in the case of polymer 
dynamics! 3 ^. With the aid of the expansion of <t+, all of models are reduced to 

d T r(z, r) - 8 3 z r(z, r) = e[d 2 z - df]r(z, r) + £(z, r) (4-3) 

in weakly stable region. When we start from (2-7) the scaled variables are given by 

r = e 3 /3t, z = P~e(x + c t), £(z, r) = e 3 (3f n (t) with e = ^^v 73 ^ an d P ~ ' 



3c ' C 

The solution of (|4-3j) is given by 



r fc (r)~/ dsexp[\ k (T - s)]£ k (s), (4-4) 
Jo 



where A^ = ik 3 — ek 2 (l + k 2 ). Thus, we obtain the correlation 

< f*(r)f_*(0) >= 2effe2( i + fc2) ex P [A fc r], (4-5) 

where 2 is the system size in this unit, and we use < Cfc( r )Cp( r/ ) >= — ^fc+p,o^( r — r ')- 
Substituting (|4-5[) into ^>(z, z', i) =< (r(z, t) —r(z', t)) 2 > which is the continuous 



limit of the scaled 4> nm , we obtain 

= 2£> G * + ^ E fc2(1 1 +A;2) {le*** - e^'| 2 + 2(1 - e M) e ^-*')} . 

(4-6) 

With the aid of , „, — -ttt ~ / and cos (^( z — ~ i l\z — z'\, the 

^ fc 2 (l + A;2) 3 ^ fc 2 (l + & 2 ) 3 I I. 

first term in the summation in eq.fl4-6|) is reduced to 



I ikz ikz' I 2 



On the other hand, the second term in the summation in eq.( |4-6| ) becomes 



where we use the approximation A^ ~ i/c 3 . Replacing the summation XmLi to the 





integral dk, and substituting (f4-6|)-(4-8) into (4-2), we obtain 



rl nt.2 r)t,2 

5 fc (r) ~ 2 / dwexp[-D G A: 2 T - —— w T l / 3 h{u)], (4-9) 

Jo 2e 7re 

where = \z — z'\, u = xt -1 / 3 , and the argument of S k is replaced by the scaled 
time. Dq is the diffusion constant for the center of mass in (4-3), and h{u) is 

,/ s f°° cos(Qu) , sin(Qu) sin(Q 3 ) . . . 
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where Q 3 = /c 3 t, and w = \z — z'\. 

In the long time limit, eq.(4^9) is governed by the diffusion of the center of mass. 
For wide range of time, however, the contribution from the second and the third 
terms in (}4-9|) are dominant because of 1/e ^> 1. In such the case, the first term is 



negligible, and h{u) converges to 

m= i dQ —^ LL= h ix i^ = wm' (4 - n) 

as time goes on. From lim^oo dw exp[— Dk 2 w/(2e)] = 2e/Dk 2 , we obtain 

Dk 2 

in intermediate time range. In the limit of small r, Sk(r) oc 1 r(l/3) ^ 2t1 ^ "■ 

Thus its Fourier transform, which is nothing but the power spectrum Pk(f) =< 
\pk(f)\ 2 > obeys 

Pk(f)~r a , a = 4/3 (as / - oo), (4-13) 

where use was made of f^° 0O dTe l27r f'[]T\ 1 / s oc /~ 4 / 3 . The value 4/3 is identical to 
the one obtained by the experiment^-* and numerical simulations Thus our 
model (2-7) r^mduces_Q; = 4/3. This result should be valid even when we start from 
fluid modelsH3)'EJ)'t3)'ll3 since the result is determined by the universal feature near 
the neutral curve as in Q4-3|). It should be noted that the appearance of this power- 
law form in the original model (2-7) is only for f < C since we eliminate the fast 
decaying mode <r_ in our analysis^ This tendency is also observed as the higher- 
frequency cutoff in the experiment H-iP . Thus, / _4//3 law is determined by short time 
behavior of the dynamics of density waves induced by the noise, which is essentially 
determined by the linear dispersion relation ~ ik 3 . 

In this section we have confirmed the universal law P(f) ~ / _4//3 of in the fre- 
quency spectrum of density correlation function from both the simulation and the 
theory. We have also clarified the mechanism to emerge /~ 4 ^ 3 spectrum which is 
related to the critical slowing down of the jiensity fluctuations. It should be noticed 
that the continuous increase of a in from a = to 2 with the particle den- 

sity is consistent with 4/3 law and our picture, because the spectrum determined 
by the noise in linearly stable uniform state far from the neutral curve should be 
white (a = 0) and the effective exponent of the power-law becomes large when the 
exponential decay (i.e. a = 2) in the off-critical region exists. There is, however, 
discrepancy between our results with the one on the experiment in liquids E^. The 
reason of this difference should be clarified in the future. 



§5. Concluding Remarks 



As we have seen in section 3, our theoretical analysis gives very precise results 
on the phase separation between jam and non-jam phases for pure one-dimensional 
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models. Of course, we do not think that our analysis is perfect. Since, for example, 
we omit the time evolution of reduced dynamical models, we cannot explain the rea- 
son why the linearly unstable branch of the coexistence curve is stable in simulation 
(see Fig.l). The validity of choices of ( |3-18| ) and ( p-23j ) are also not confirmed from 
mathematical point of views, although these choices work very well. To clarify the 
above points will be future subject of the research. 

Let us comment on the universality class of traffic flows and granular flows. 
All of models introduced here except for ( |2-3| ) have asymmetric kink-antikink pairs 
and qualitatively resemble behaviors with each other. In fact, wet^P have already 
checked the Quantitative validity of our methods presented here for the fluid model in 
traffic flows Eif. The results are almost identical to those explained in section 3. On 
the other hand, OV model in fl2-3| ) which is a special case of the above generalized 
models loses some universal properties. It should be noticed that the limitation of OV 
model has been suggested by Komatsu and Sasatf (see the last part in their paper). 
Therefore, we believe that our analysis is meaningful to characterize universal feature 
of one-dimensional dissipative flows such as granular flows and traffic flows. 

On the other hand, we have extended a one-dimensional deterministic model to 
a stochastic model supplemented by the white noise. This model clearly reproduces 
P(f) ~ / _4//3 l aw as i n the experimentliiP. We also give a simple argument for 
the reason why we obtain / _4//3 law. Readers may be skeptical whether our one- 
dimensional models supplemented by the white noise can describe true behavior 
of quasi one-dimensional systems in spite of good agreement with the experiment. 
To reply such the question, we have already introduced a simple model for two- 
lanes traffic flows©, where the vehicles have spin variable to specify what lane is 
chosen and lane-change of vehicles is governed by Glauber dynamics for anti-ferro- 
magnetism. Our preliminary result of simulation also supports /~ 4 ^ 3 law. Thus, 
/ _4 / 3 law is believed to be universal without regard to the choice of models. The 
effect of nonlinearity in the theoretical argument and the validity of the introduction 
of the Gaussian white noise in (2-7) will be discussed elsewhere. 

In conclusion, we have proposed a simple generalized optimal velocity model 



(2-1). Based on the perturbation analysis of asymmetric kink solution (3T7) we 
obtain the selected values of its amplitude, propagating velocity and width of kink as 
in ( |3-34| ) . The accuracy and relevancy of the solution has been confirmed by the direct 
simulation. We also extend a one-dimensional model to a quasi one-dimensional 
model with adding the white noise. Prom the simulation of the noise sustained 
model and analytic linear relaxation theory we have confirmed the universality of 
/" 4 / 3 law in power spectra in density auto-correlation function. 
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